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In many instances one 3ncoT;!ntc'''s 'vitb the problem 
of flow through non-horaogeneous and anisotropic porous 
media. In general, for analysing such flow situations 
either numerical or analog methods are used. The present 
study deals with the analog solutions for the problem of 
confined flow below the weir resting on non-homogeneous 
and an i so tro pic soils. ' " 

The non-homogen ity considered in the present analysis 
corresponds to the jointed rock mass which is characterised 
by two permeabilities, namely primary and secondary per- 
meabilities corresponding to rock material and joints, 
respectively. ,v.-'$he analysis is presented for a weir with a 
vertical sheet pile resting on an infinite rock mass consider- 
ing one major joint either to be vertical or horizontal. The 
influence of (i) disposition of sheet pile (ii) length of 
sheet pile, (iii). depth at which horizontal joint is 
located &(iv) the presence and location of vertical joint, 
on the equipotential lines and the uplift pressure is invest- 
igated. The results clearly bring out the necessity of 
0 )nsidering the disposition of joint in the analysis of such 
a flov\f situation.v ■ 
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Similar to nor-homogenity , the anisotropic nature 
of the soil also governs the confined flow behaviour 
considerably. In general, by transformation technique, 

-the anisotropic flow domain can be converted into an 
equivalent isotropic region for determination of uplift 
pressure, exit gradient etc;> ^^alking the general case 
of anisotropy and varying the lower flow boundary condition, 
analog solutions have been obtained for a weir with one 
sheet pile.^' The influence of (i) the location of sheet 
pile (ii) inclination of major axis of ellipse of ani- 
sotropy with respect to horizontal direction, (iii) lower 
flow boundary conditions (either highly pervious or 
impervious) is investigated in detail, ■y'' 

As it is very difficult to obtain the solutions for 
a two layered soil system which is not parallel to 
horizontal direction, it is hoped that the results pertain- 
ing to the different lovrer flow boundary conditions would 
be useful as lower and upper bounds for a layered soil 
system seperated by the geometry corresponding to the 
flow boundary. 
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INTRODUCTION 

Seepage oompui;ations play an important role in. the 
design of hydraulic structures, retaining structures and 
stability of slopes* In case of design of hydraulic 
structTures the specific problems to be dealt with can be 
classified as: 

1. Estimation of q.uantity of seepage and zone 
of seepage. 

2. Uplift pressure distribution at the base of ihe 
structures, and 

3. Estimation of exit gradient to safeguard 
against piping. 

It should be noted that in many instances the 
damage occurring in hydraulic structures can mainly 
be attributed to the destructive effect of seepage. 

Even though the basic concepts of theory of fluid 
flow through porous media have been expounded nearly a 
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century back, it is only since th.e beginning oi 
twentieth century the subject has received a scientific 
treatment. of the reasons for chgq.uered development 

of the subject in the early stages has been the hetero- 
geneous nature of the soil at microscopic level. 

However, now it is well^accepted fact that most of the 
seepage and groundwater flow problems can be analysed 
on the basis of larcey's simple linear law proposed in 
1856. In addition to the proof of validity of Darcy's 
law, the developments in hydraulics and hydromechanics 
have made it possible to derive several theoretical 
solutions which are applicable to many field situations. 

Most of the available theoretical solutions of 
groundwater flow assume tho porous media to be isotropic 
and homogeneous with respect to coefficient of permeability* 
However, materials occurring in nature are seldom 
isotropic and even less frequently homogeneous. Ihile 
many adequate general methods of solving seepage problems 
through isotropic and homogeneous porous media have 
been developed, problems of confined flow under weirs 
resting on anisotropic and non-homogeneous soils 
have received little attention due to associated ■ 
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mathematical difficulties. Plow through, anisotropic 
porous media is generally analysed by first transform- 
ing the anisotropic actual flow domain to a fictitious 
isotropic flow region by a suitable co-ordinate trans- 
formation and then solving the probloms for the 
equivalent isotropic system. Por the problem of confined 
flow below hydraulic structures resting on non-homogeneous 
soils, theoretical solutions are only available for the 
casB of a two layered soils system wiih. both the 
layers having same thickness. However, in nature the 
non-homogeneity can occur in the form of localized 
pockets of impermeable or permeable soil and ^jointed 
rock forms having primary and secondary permeabilities. 

Por the foremen tioned situations no theoretical 
solutions are available in literature because of the 
complex flow boundary conditions at the interface 
joints. 

Review of the availa.ble literature neveads that 
several methods have been used for finding solutions 
to problems of flow through porous media. Some of the 
commonly adopted methods are! conformal mapping 
finite difference method, finite element method, graphical 
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approach, energy methods and analog methods, ilany 
probleras vdaich are difficult to be solved analytically 

can be relatively easily solved by electrical analog 
method. It has been found that the accuracy of the. 
results obtained by analog method is good enough for 
many engineering problems. 

As has been stated earlier, the available knowledge 
regarding the confined flow problems for hydraulic 
structures resting on anisotropic and non-homogeneous 
soil is far from complete. With this in view the present 
investigation is primarily concerned with the study of 
some cases of two dimensional confined flow under weirs 
resting on non-homogeneous and anisotropic media, using 
electrical analog approach. Also the present study 
primarily confines only to the distribution of uplift 
pressure at the base of h,ydra.ulic structures. 

The scheme of presentation in the thesis is as 
followss- 

A review of the literature in the field of seepage 
in anisotropic and non-homogeneous porous media and use 
of electrical analogy relevant to the present study has 
been presented in Chapter 2. 
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Chapter 3 deals with the elements of electrical 
analogy, its usefulness and its limitations for solving 
problem in groundwater and seepage. Details of the 
analog field plotter used in the prestint investigation 
are also presented. 

Chapter 4 deals with the problems of confined 
flow through the jointed rock mass which is a particular 
type of non-hoinogenasiity commonly encountered. 

Dor this case it is difficult to arrive at exact 
solution due to existence of two permeabilities, the 
primary permeability of the rock mass and secondary 
permeability of flow throu^ joints. Recently Madhav 
and Lakshmidhar (1968) presented analog solution for 
flow below the weir founded on a jointed rock mass, lo 
study the effect of jointing the joint was either 
considered to be horizontal or vertical and results thus 
obtained qua.litatively brought out the effect of such 
major joints on the flow behaviour. In the present 
work analysis is presented for a weir with a vertical 
sheet pile fDunded on jointed rock mass. Dor simplicity 
the joint is assumed either vertical or horizontal (Dig. 1(a) 
and 1(b)). As the effect of joint has been simulated by 
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drawing a strip of silver paint on analog paper, the 
results are once again qualitative than quantitative. 
However, in the absence of any available solutions the 
analysis will serve to bring out the influence of joint 
position as affecting the uplift pressure. Results are 
presented in non-dimensional form for various locations of 
joints, length and disposition of sheet pile. 

The effect of flow boundary as affecting the uplift 
pressure for a weir resting on anisotropic soil is 
analysed in Chapter 5* In iiie initial portion of the 
chapter, review of the transformation technique is presented 
which transforms the anisotropic region into an equivalent 
isotropic region. In the later part of the Chapter 
results of the study bringing out the effect of lower flow 
boundary are presented. Various cases analysed are: 

(i) Confined flow under a weir on anisotropic 

soil of finite depth with sloping impervious 
boundary (ff’ig. 2). 

(ii) Confined flow under a weir resting on an 
anisotropic soil of infinite depth with 
sloping pervious boundary (Rig. 3). 
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(iii) Confined flow under a weir on an 

anisotropic soil of finite depth under- 
lain by a draining layer and with sloping 
pervious boundary (lig. 4). 

Por all the cases considered, the effect of various 
parameters as affecting the uplift pressure is presented 
in the form of non-dimensional charts. It should be 
noted that for various cases considered, no theoretical 
solutions are available at present and as such the 
analysis would be useful while dealing with problems 
encountered in field with such complicated flov; boundary 
oonditions. Also as it is extremely difficult to 
obtain theoretical solutions for a layered soil system 
which is not parallel, the solutions can be thought 
of either lower or upper bound to a two layer soil system 
separated by the flow boundary which is either assumed 
pervious or impervious. 

The last Chapter (6)dGals with the general 
conclusions based on tie present study and also includes 
suggestions for further Research. 
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CHAPTER 2 


EBTISW OP LITBItiiTUHB 

2.1 GEiraRAD; 

In this Chapter a brief review of the literature 
pertaining to the two-dimensional confined flow problems 
is presented. In the first section, review is presented 
of the various analytical investigations dealing with 
the problem of weirs on permeable foundations. More 
details are given of the recent studies which include 
the effect of anisotropy and non-homogeneity in the 
analysis. In the second section of the chapter, review 
of relevant electrical analog studies is presented. 

2.2 THE ORE TIC AI SOLu'TlOhS: 

Prom the time Pa:^'A( 16356 ) established the linear 
relationship between the seepage velocity and hydraulic 
gradient, many theoretical solutions have been presented 
for both confined and unconfinad problems. Some of the 
notable contributions which have a direct bearing on the 
present investigations are: the case of confined flow 
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under a weir with, a vertical sheet pile and resting 
on porous media of infinite depth (Khosla, etal 1954), 
confined flow under a weir with a. vertical sheet pile 
and resting on a finite porous layer underlain by 
either a impervious stratum (Muskatf 1937) , Pavlovsky 
( 1933 )) or a horizontal draining layer (G-rinsky 1950) , 
confined flow under a depressed weir (Pavlovsky 1933) 
the case of an inclined sheet pile in a semi-infinite 
horizontal porous medium (Vergin, 1940) and the case 
of a confined flow under a weir on an anisotropic medium 
of finite and infinite depth (Reddy and Mishra 1972). 
Theoretical solutions co.nsidering the non-homogeneity 
of the porous medium are presented by Polubarinova 
Kochina( 1962 ) for the cases of vertical sheet pile and 
a weir without sheet pile, resting on, a two-layered soil 
having individual layers of equal thickness. 

In all the above cited works solutions have been 
obtained by me.ans of Schwartz Chrtstoffel transformation 
technique. The final solutions are obuained in general 
in terms of hyperbolic functions, or elliptic functions 
or incomplete and complete beta and gamma functions. 

Most of the above theoretioal solutions are based 
on the assumption of isotropic and homogeneous soil medium. 
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However, in reality, one often encounters situation 
where the soil is either anisotropic or non-homogeneous 
or both- 

A general the'^retioal analysis of flow of fluids 
through anisotropic soils has boon presented by Eerrandon 
( 1948), Scheidegger (I957j, lolubarinova-Kochina (1962), 
Harr ( 1962), I/l&orcus ( 1962) fund more recently by Misra 
(1972). 

In general for taking into account the anisotropic 
nature of the soil, the actual flow do main is transformed 
into an equivalent isotropic flow domain by co-ordinate 
transformation. The details pf transformation can be 
found in the book by Harr (1962), Polbarinova - Kochina 
( 1962), and le^wiest ( 1965 ) . By mcking use of ttie fore- 
mentioned transformation technique, Virgin (1940) has 
solved the problem of a vertical sheet pile in an 
anisotropic medium. Adopting essentiailly the same 
approach (Reddy and Misra 1972) have solved the problem 
of confinod flow under a weir on an anisotropic medium 
of finite and infinite depth. 

Also in many inscances, the region in which the 
seepage takes place is essentially of a non-homogenoous 
character. In such conditions seepage problems become 
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complicated. For example problems of seepage beneath 
water retaining structures often require the consider- 
ations of layers of soil with different permeabilities 
which also may be anisotropic. The permeability being a 
function of void ratio, may decrease with depth within 
a layer and stratum may be lense like. If the non-homoge- 
neity. can be characterised by stratification consisting of 
sequence of thin homogeneous isotropic layer of different 
permeability and thickness, the problem can ultimately be 
converted into a form that is tractable by one of the 
technique listed in Chapter 1. In general this is not the 
case and there are fewer methods available for the solutions 
of irreducible non-homogeneous, anisotropic seepage 
problems. For simple flow situations Pavlovski and 
Kamenski have given some approxim'.. te methods for reducing 
a multilayered soil to equivalent homogeneous one. 

'Flow net skotciiing can be done in practice, but 
it is not practicable except where trie non-homogeneity 
is relatively simple. The Kele-shaw model can be used 
for non-homogeneous isotropic conditions but requires 
carefully constructed appratus- The most suitable methods 
in current use for solving more general problems are 
electrical analog and numerical methods. The electrical 
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resistance network Sc otH: ( 1963 ) , Herbert and Huston 
(i 960 ) provides a tecnriique of considerable versatility 
whicli may also be used for 3 - dimensional problems. A wide 
variety of problems can also be solved numerically using 
digital computers. Ihe finite elements method 
2 ienkiewic 2 ( 1966 ) can be effectively used for analysing 
the problems of flow through anisotropic and non-homog- 
eneous soils. 


Flow through jointed rock is also a case of flov; 
through non-homogeneous porous media due to occurrence 
of two permeabilities of rock mass. Kransnopol' skil 
( 1912 ) proposed a formula for seepage in jointed rocks, 
according to which the head loss is proportional to 
the square of seepage velocity. 


A similar formula was derived by Puzyrevskii 
( 1930 ) for soopage in a stone talyis. 

Serafim and Oampo (1965) calculated the average 
permeability Keiv for the rock mass as: 



_AiL ( £ 

12 V i 



d 


width of joint; 
spacing of joint; 
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"Yuj = unit weight of fluid; 
and = viscosity of the fluid. 

The above equation is useful when one is interested 
in calculating the quantity of seepage through the rockmass, 
hut estimation of pressure head at different location is 
not possible with the above equations. 

¥ittke and Louis (1966) neglect the primary 
permeability of the rock mass and base their analysis 
on the geometry of the joints. Flow equations are develop- 
ed by Louis (1966) for seepage through open or filled 
joints considering the form and the surface roughness 
of the joints. 

For quantitative study finite element methods 
(Zienkiewicz ( 1966) )ha.ve been more readily applied to 
the seepage through rock foundation. A simple numerical 
procedure has been recently developed by Madhav (19?©) to 
make quantitative study of flow through jointed rock mass. 


2.3 ELSCTEECAL iuvjJjOG- SOLUTIONS: 

From analogy between Ohm’s Law for condition 
of electricity and Darcy’s Law lor flow of water through 
porous medium; Pavlovsky suggested as early as in 1921, 
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that an analog method oould he employed for the purpose 
of studying the uplift pressure under dam. Some 
attempts were made in this direction hy 6c Coworkers 

( 1954 ) in America hut no definite conclusion were arrived 
at. However, suhaeg.uently the analog method was improved 
and was successfully used by many research workers notable 
amongst them being (1935), Wycoff and Seed ( 1935), 

Selim ( 1947 ) and Vaidhianathan ( 1955). 

For the first time three dimensional electrical 
analog technique was used hy Reltov in 1936 to study 
the flow pattern under a hydraulic structure built on a 
heterogeneous soil mass. 

H^r-t ( 1935 ) employed this technique to investigate 
(i) uplift pressure (ii) exit gradient and (iii) quantity 
of seepage under dam on sand. Electrolytic tank of size 
24 X 26" with conducting media as salt solution was 
used for conducting experiments. Also the results of 
his study showed remarkable closeness to the already avail- 
able theoretical solution, thus proving the usefulness of 
analog technique. 

Eor the analysis and design of weirs resting on 
permeable foundations ¥aidhinathan et al (1955) used 
extensively the electrical analogy method. They studied 
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•the uplift pressure distrilutiDU for liie cases: 

(i) A simple impervious floor flush with the. surface of 
the pcrous strata. 

(ii) The same case with sheet pile. 

For both cases, theoretical solutions and direct 
measurement in hydraulic models existed, so the results of "the 
electrical analogy method could be compared with direct 
observation of pressure and with theoretical solutions. ■ 

In all the cases the agreement between electrical analog 
solution, theoretical solutions and observed direct measure- 
ments of hydraulic models was reasonably good. Subsequently 
Vaidhianathan (1955) presented results of exhaustive 
investigation pertaining to "the weirs with number of sheet 
piles resting on finite and infinite porous stratum. In 
addition to the above mentioned problems Vaidhinathan (1955) 
successfully used this method - for solving problems relating 
to depressed weirs, tube wells, coffer dams and earthworks. 
The appratus employed by Yaidhi wreath an, for the study 
consisted of two parts the nodal tank (4'x3' x 3”) and 
the source of alternating carrent. Two conducting plates 
of copper each 1.5' long respresenting the surface of the 
porous strata on the upstream and downstream of -the dam 
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were used. The base of the dam was simulated by 

1' long elJOnite plate, flush with conductors. The 
wires 

potentiometer / served as a potential . divider, 
luthra and Joglekar (1961) employed electrical analog 
technique to investigate the distribution of uplift 
pressure below hydraulic structure for the situation 
when the sub-soil thougn. permeable, is in stratified 
layers. They considered the following aspects of 
stratification: - 

(a) Oases where a number of layers of different 
permeability exist in subsoil. 

(b) Oases where an impervious stratum exist under i 

a pervious sub-soil medium. 

(c) Oases of an..impervio us layer of irregular shape 
and localised pockets of impervious soil. 

Meehan, R.L. and Morgenstern, N.R. (1968) employed th 
technique for getting approximate solutions of seepage 
problems for following cases 

(i) Sheet pile cutoff in homogeneous isotropic soil. 

(ii) Sheet pile cutoff in layered system, and 

(iii) Earthdam on anisotropic layered foundation. 



17 


They employed a conducting resistance networks 
which was prep.ared hy drawing lines of graphite 

ink on a. sheet which conducts electricity when dried. 
Permeability was varied h'j increasing or decreasing tbe 
spacing of lines. 

Stefen, H, and Ooworker (1964, 1970) solved the 
problem of two dimensional seepage under a flat-bottomed 
structure with a sheet pile cutoff Vv'all at upstream end 
resting on two horizontal layers of different permeabilities 
and of eq.ual thickness. Conducting paper was used 
to represent flow field. Since the permeability is 
proportional to the number of layers, it is possible to 
adjust permeability locally by using several layers. In 
order to vary permeability ratio which ranged from 
0.25 to 8.00, several sheets of paper were glued together 
with amylaoetate. The experimental technique was 
verified against the analytical solution for seepage 
flow under impervious structure without a cutoff wall, 
as given by Harr (1962). 

Per the analysis of confined flow below a flat 
bottoned weir resting on jointed rockmass Madhav and 
lakshmidhar ( 1968) used the electrical analog method. 
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Analog field plotter was employed in the study. 
Conducting graphite (Teledeltos) paper was used for 
drawing equipo tential and flow lines. Primary permea- 
bility of the rock mass v/as simulated by electrically 
conducting graphite ( L-eledeltos) paper while secondary 
permeability by a thick line of silver paint on the 
paper . 

Raza ( 1971 ) carried out investigation to 
evaluate the distribution of uplift pressure below a 
dam and to determine the effect of upstream sheet piling 
and presence of localised non-homogeneities on the 
uplift pressure. He also studied the effect of size, 
shape and number of discon tinutie s present in the region 
below a dam as affecting the uplift pressure. 

Punmia cndKhullar ( 1972) presented results of 
analog study dealing with uplift pressure below apron 
founded on pervious medium of finite depth with 
downstream cutoff under scour conditions. 

It can be seen from the above brief review 
that tho electrical analog method is one of the 
varsatile techniques which can be employed for solving 
problems in the area of seepage and ground water. "With 
this in view in the present investigation electrical 
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analog tochniq.U0 is used to solve some seepage 
problems talcing into account anisotropy and non-? 
lioraogeneity for which neither analytical nor experi- 
mental solutions aro available so far- 
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EIECTEIOAL ANALOGY 


3.1 GBEBEALj 

In many branches of engineering and physics, 
Laplace's equation is encountered as the field 
equation. Among these areas are ; the steady state 
flow/ of groundwater, electricity and heat. The basis 
of analogy becomes obvious when one considers the 
governing laws in various areas, that isj Darcy's ].cw 
of steady state flow; Ohm's law for current 
conduction and Fourier's law for heat conduction. 
Theoretically, any of the above analogie s can be 
used for analysing the problem of confined flow. 
However, building up of a thermal analog model to 
simulate soil Wi/ater flow situation is quite a formid- 
able task since it is difficult to form boundaries 
under the requirement of ihermal insulation along 
the flow boundary' and the maintainance of constant 
temperature along the potential boundaries (Scott, 
1963). Because of this, thermal analogy is not used in 
practice. However, one can conveniently use 
electrical analogy to solve the problem of steady 
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state flow through porous medium. Also Hie 
electrical analog has added advantage of simplicity 
of experimentation and the solutions can be obtained 
with ease and rapidity. The correspondence between 
the steady state flow throuA porous media and the 


steady flow of electric current 
form below: 

Steady State Seepage 
Darcy’s law for water seepage 
Q = (KA. h/l) 

1. Q = quantity of discharge 

2. K = Coefficient of 

permeabili ty 

3. 1 = Gross Sectional Area 

4- h = Pressure head 

5. 1 = Length of Seepage path 


is presented in tabular 


Electric Current 

Ohm' s law for conducticrA 
of electricity 

I = (CA' V/l' ) 

I = rate of flow of 
electricity 

C = Current Conductivity 
Coefficient. 

A’ = Cross-sectional 
Area 

Y = Electrical 

Potential (Yoltage) 

l’= length of path of 
electric current 



22 


Field Equations: 

6.. 9^ ii = 0 V^y = 0 


7 . 


Impervious 

B h 
■aiT 


loound ary 
= 0 


Insulated boundary 

M =0 


In case of heat flow through conductors, or 

that of electricity, an analytical treatment is 

possible. This is because the conductors in these two 

cases are generally metals or crystals and they have 

regular isotropic or anisotropic structure. In contrast, 

the geometrical boundaries surrounding the intersticial 

water in sand or caly, always differ and do not repeat 

in any, definite arrangement. Also, the nature of 

packin g changes the transmission i.e. changes the 

conductivity. Thus, tho analogy, therefore fails in 

the iegions of . molecular or microscopic behaviour 

porous 

of flow of water through/media . It can only be viewed 
as applicable to macroscopic flow ( Vaidhinathan 1955). 
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3.2 EIECTEICiL MUJjOG- MOHELS ; 

For solving the problems of confined flow 
through porous media both continuous and lumped 
parameter amalog are employed. The details of some 
of the commonly used techniques are given belows 

( 1 ) Gon tinuous Bfodels; 

(a) Electrolyte . Trough Me thod i 

This has been used quite extensively to solve 
two or three dimensional seepage problems. The 
electrolyte or trough consists of a layer of current 
conducting electrolyte contained within a suitable 
insulating trough whoso size and shape are regulated 
to conform with the outlines or boiindaries of the 
problem to bo studied. Current is fed into the 
electrolyte through suitably immersed eloctrodo, and 
oquipotentials linos are traced with a detecting 
instrument. A schematic diagram of a typical electro- 
lyte trough model is shown in Pig. ^.1. 

(b) Analog Field Plotter - 

m 

It consists of a thin shoet of electrical 
conducting paper (0.004" thick) in which an electrical 
current flow pattern is set up by means of suitably 
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attached aaad energized electrodes. The resultant 
potential-drop established by this current pattern is 
detected and marked directly on the paper by means of a 
searching sytlus in conjunction with a high sensitivity 
detecting instrument. Since the nlow of current is in 
two dimension hence this metbod is only applicable 
to the analysis of two dimensional seepage problems. 

(2) Lumped Parameter Model ;- Fig. ^.2(a) represents 
a section of soil medium throu^ which steady state . 
seepage is taking place. The actual flow domain has 
been discretized in finite number of points as shown in 
Fig. 2(b) for developing an eq.uivalent analog. As can 
be seen from Fig. ^.2(h) in the corresponding analog 
the continuous nature of the medium is approximated 
by ai, number of discrete resistances joining the points 
a, b, c . . . etc.j corresponding to A, B, C . . etc. 

in actual flow domain. Thus the continuous conductance 
AB, BC, . . . etc. are replaced by lumped resistars ab, 
be . . . etc. , the analog electrical conductors being 
related to the hydraulic conductance by a fixed scale 
factor. Thus, in lumped parameter model a permanent 
board set with connectors, usually in square array 
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at a suitable spacing is used to solve the confined 
flow problems. Detailed account of luijped parameter 
analog can be found in the book by Scott (I963). 

In comparison to analog field plotter 
electrolyte trough method has the disadvantage that 
it requires rather elaborate and laborious initial 
setup. In addition the tochnique of conducting 
experiments may involve an appreciable amount of 
polarization and phase shift in certain cases, and has 
the inherent difficuliyof using a spillable electrolyte 
whose electrical properties are sometimes difficult 
to regulate and maintain uniformly constant over a 
period of time. 

In the case of lumped parameter analog model, 
the network can be extended without any difficulty 
to radial, two dimensional or even three dimensional 
problems. However, in case of continuous analog 
employing dry electrically conducting apaper the 
analysis is possible only for two dimensional 
problem. Also inclusion of non-homogeneity is relvctively 
easily accomplished in case of lumped parameter analog 
when compared to continuous analog. Transient 



hydraulic flow can also be simulated in a lumped 
parameter analog by connecting an electrical capacitor 
to each node point of 'tlic ;inalogous electrical net. 

In the present investigo'^? on rs the problems 
of confined flow to be investigated could be easily 
analysed on a continuous model, using dry electrically 
conducting paper, the seme has been used for various 
studies. The details of analog field plotter used are 
given in following section. 

3.3 MIOiOQ FIELD PLOTOERi 

A schematic diagram of the tmalog plotter is 
shown in Fig. 2.3* Plate 2.1 shows tho photograph 
of the complete setup. The various components of the 
an alog plo tte r are 5 

( A ) Conducting Paper Sheet; 

The sheet of conducting paper in which the 
electric current antilogy field pattern is established 
consists, of a piece of type 'L’ teledeltos paper 
apprcximately 0,004” thick. The conducting paper 
having resistance of approximately 4000 -^^/sq,. has 
been found very well suited to the use in the field 
plotter. The ratio of the paper resistance ( 4 OOO "'^/sq, ) 
to silver painted electrodes areas (1 to 4~’h/sq) is 
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adequately high to permit the conducting of the 
element area to be considered a.s being infinite 
in most of cases. The average peper resistance uni- 
formity of +3 ¥ to + is believed to be 

adequate for hy-f-eir the majority of field plotting 
problems. This nonuniformity property in resistance 
appears to be a permanent property introduced in the 
processing of paper. 


Base Board; 


The base board used as a plotting base may 


consist of an ordinary 23” x 31" drawing board of 

3 " 


4 


thick board of insulating material, preferably of 


plywood, which has a smooth flat surface on which the 
conducting sheet is mounted. 


( ^ ) I nstrument Null Detector: • 

It serves to indicate the location of the 
poltite of zero potential which serves to mark a 
particular eqnipo tentiol line being traced. The 
voltage level of the particular equi -potential 
being traced, is determined by the particular slider 
position chosen on the voltage-dividing potentiometer, 
the voltage level being adjustable in small 
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increaments anywhere he tween 0 ond 100 of the voltage 
applied to the pl'^ttor. 

( ^ ) S g arching; and Markina; Stylus; 

The se 'inching probe serves both to explore and 
locate tho zero potential halanco points comprising 
the equipotential lines and as a marking stylus whose 
point is used to make small indentations or perforation 
in plotting surface for recording the equipotentials . 

The stylus consists of a moderately sharp pointed 
brass rod with an insulating sleeve, about 5” long. 

It is connected to the terminal marked P, which is the 
slider circuit of the instrument voltage divider. 


(E ) Voltage. Divider; 

Tho type of voltage-divider unit multiturn 
potentiometer has been found to be very satisfactory for 
use with field plotter, having the marked advantage of 
high precision coupled ?ath very fine incremental 
adjustment of the slider. 

( ^ ) Electrode and Boundftjy Control MaterigJ- -& Method 
of Poruiinjg: Silver Painted Electrode; 

Among the various miscellaneous and supplies 
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wliich are used at various times with the field plotter in 
setting up problems are the followings 

(1) Silver conducting paint 

(2) xicetono ttij-uol for reducing silver paint to 
proper consistency. 

Since it is essential that the electrode areas 
be generally of very hi^ conductivity, copper wires 
have been attached to the left and right-hand electrode 
strips inorder to insure that these strips are lov/ 
enough in resistance to represent infinite equivalent 
conductivity as far as plotting circuit is concerned. 

Method of foiming silver paint - Low resistance 
electrode areas may be conveniently established on 
conducting paper sheet by moans of silver conducting 
paint. Silver conducting paste, which comes with about 
the consistency of molasses, is thinned from 1 to 3 
Piirts by volume of a thinner solvent such as acetone. 

In special cases where exceptionally low- 
resistance electrode areas may be required a second or third 
coating of silver paint may be applied after the first 
coat has dried or both sides of paper may be painted 
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ovor the desired electrode areas. It usually takes 
24 hrs to get the paint dried. 
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CHAPOEE 4 


GOWPINEP PLOW miER A WIR EBSTING ON A o^OIFIED ROOK MASS 

4.1 INTROPUGPION: 

In many instances liydraulic ctructuree ■^"^unded on 
rock-mass which hy definition consists of rock material 
and planes of weakness due to geological separations. For 
analysing the flow under such hydraulic structuiss, it is 
essential to know characteristics of fluid flow throu^ 
jointed rock mass. It is obvious that because of the 
presence of joints, either major or minor, the flow domain 
can no more be idealised as homogeneous. Ono of the ways 
in which the non-homogeneity can be ta'-^en into account is 
hy characterising primary permeability for rock material and 
secondary permeability for joints. By its inherent nature 
the former is small while the secor>dwf ’figTO ability by 
comparison a much larger parameter. Based on this approach, 
recently Serafin and Gampo (-1965) and 1^?lttko and Louis (1966) 
have presented basic approaches to the analyses of flow 
through jointed rock mass. Recently in 1968, Madhav and 
Lakshmidhar developed electrical analog method to study 
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(i) primary permeability for the rock material and 

(ii) Secondary permeability for joints are considered in ^ 
the analysis. As the end result of any such analysis 
would considerably depend on the values of these permeabi- 
lities, in the following paragraphs further details 
regarding these parameters are presented. 

The primary permeability of rock material bounded 
by either major or minor joints is due to the presence of 
in rock matrix. Por finding the primary permeability 
in the' laboratory, experiments are conducted on cores from 
rock material either using gas or water as the fluid, 

(Stagg and Zienkewicz). In general the primary permeabi- 
lity of rock mass is of the order of 10 to 10 ^ cm/sec. 
(Serafin and Oampo, 1965). Also the primary permeability 
is anisotropic in nature being more in horizontal direction 
compared to the vertical direction. 

Terzaghi ( 1962) defines the secondary permeability 
as that resulting from the flow through open smd conti- 
nuous joints in the rockmass. It is obvious that by 
very definition the secondary permeability of a rock 
mass ’Will be dependent on the opening and spacing of 
joints. In case of clay filled joints, tbe permeability 
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will be a. reduced one than the distance between two 
joint surfaces. For a. rock mass with a system of 
parallel joints of width e, separated by a distance, d, 
the eq.uiTalent permeability (K) of a porous body is 
given by Serafin and Compo ( 1965) as 


K 


F 


^ w e 
12 1-i 


(4.1) 


where I' and Y:*j are coefficient of viscosity and 
unit weight of the water. Since the above eq.uation 
is independent of the geometry of the joint pattern, 
it would be useful only for c.alculating th-e quantity 
of seepage through rook mass and not for estimating 
the pressure head at different locations. 


The coefficient of secondary permeability 
can be determined by a "Joint permeability" test and 
expression for secondary permeability (k;s) can be 
written as: (Serafin ?;ndOampo, 19 65) 


Kp = 



(4.2) 
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Even for small values of e, the value of Ks worksout 
to be mucii larger than the primary permeability (K.p). 
Depending on the type of rook and the extent of jointing, 
the ratio of K/kp can be of the order of 10 to 100 
or even more. 

Erom the above discussion, it can be concluded that 
the 2Xickmass can be idealised as a medium with a smaller 
primary permeability but divided by sets of joints 
having a high secondary permeability. For the study of 
confined seepage through such a medium with a view to 
determine the uplift pressure, very few 'Withods are 
available. A complex mathematical method has been 
developed by Barenbldtt et al (i960) which may be 
difficult for application to many field problems. Also 
the rockmass cannot be hypothetically replaced by a 
series of alternate layers havng different permeabilities 
as the joint thickness is negligible in relation to 
the rock material. Recently Madhav and Lakshmidhar 
( 1968) modified the electrical tihalog method to 
simulate the flow through jointo^d rock mass. They 
analysed the problems of a weir resting on an infinite 
rockmass with either a vertical or horizontal joint. 
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The conducting teledeltos paper served as rock material 
having privaary permeability while the lines dravm with 
silver paint served to simulc/ce secondary permeability. 

E s sen ti ally j the same approach as the above is 
used in the present study to find out the uplift pressure 
below a weir Vi/ith a vertical sheet pile resting on 
jointed rock mass. Pig. 4.1(a) and 4.1(b) show the 
schematic diagram of the problem under consideration. 

The different parameters varied during the study are 
also indicated in Pig. 4.1(a) and 4.1(b). 

Por presenting the results in non-dimensional 
form following dimensionless parameters are defined? 


/\ = 

(d/b) 

h = 

( Vb) 


(b/s) 

'VA = 

1 

(b/b) 


where 

d = distance of horizontal joint from base of 
the weir. 


b = base width of the weir. 
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1 = dist.nncG ol vorticRl joint from *0' 

(Pig. 4.1(a) in d (b). 

b^ = distanco of shoot pile from '0* 

Por various values of A 'll d and V.; the uplift pressure 
distribution at the base of tho woir was obtained so 
also the eq.uipotential lines in the flow domain. Tables 
4.1(a) and 4.1(b) summarize the values considered for 
parametoiB ,01. .|_i and A- As can be seen from tables, 
tho relevent parameters have been varied so as to cover 
m.3ny field situations. 

It should be noted that the analysis is confined 

only either to a vertical or horizontal joint. However 
the method being general enough can be easily extended to 

other configurations of single or multiple joints. The 

details of experimentation are presented in next section. 

4.3 HST'AIIS OP IlZF5RD:I55rTS: 

The analysis presented is based on the following 
as sumptions : 

(1) The joint is either vertical or horizontal. 

(2) The rock mass is homogeneous and isotroi.ic. 
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(3) Joints have been treated as tabular so that 
they can be depicted by lines in two dimensional 
represent.ation . 

The primary permeability of the rockmass is 
simulated by the teledeltos paper. To include secondary 
permeability of the jointsp lines can he drawn ohHhe 
paper with higher conductivity. This con be accomplished 
in following throe ways: 

(1) lines of uniform v^idth and thickness of silver paint 

( 2 ) Thin strips of silver foil attached to the teledel" 
paper and 

( 3 ) Thin electrically conductive wires fixed along 
the joint location. 

Out of these three methods the first one has 
been used in present investigation because it offers 
the best contact between the paper .and highly conducting 
medium. The width of the joint should be negligible as 
coiiiparod to the depth of rockmass. Since it is very 
difficult to drav^f one single line of uniform thickness 
and width, lines of 0*2 cm have been drawn. The ratio 

_3 

of joint width to depth of rockmass is 4*4 x 10 
which is q.uito small. 
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In order to reduce boundary effect', the dimension 
of the model is choosen as follows: 

(i) The ratio of the distance Ifeft either on right or 

left side to the half base width is kept 3* 5* 

(ii) The ratio of the depth of the rockmass to the 

maxinium depth of sheet pile is 6.5* 

As per latest work of Mutukuma.ratn and Kulandaiswamy 
( 1972 ), so long as the distance left on either sides 

is 2 to 3 times than half the base with of the weir and 
the ratio of the depth of rockma.ss to the depth of sheet 
pile is more than 4*0, the error introduced by boundary 
effect would be within 3 • As such, in the 

present study the error crept in due to boundary effect 
Coin be expected t'^ be within 2 io. 

4.4 RESULTS ml) DISCUSSION: 

The results for horizontal and vertical joints 
as affecting the confined flow behaviour below a weir with 
a vertical sheet pile are presented below in two sub- 
sections separately. 

( ^ ) H orizontal Joint; 

Rigs. 4*2 and 4*3 show the effect of location of 
joint as affecting the eguipo tential lines for a sheet pil- 
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at ups cream end of the weir for =0.2 and 0.6 
respectively. In both these figures for the same 
location of joint, equipo tential lines are plotted for 
various length of sheet pile. As it can be seen from 
these figures, the presence of the joint considerably 
changes the pattern of equipo tential lines. Also by 
comparing the eq.uipo ten tial lines for the same value 
of ^ in both the figures it is seen that the location 
of joint i.e. its relative distance from the base of 
the weir also influences the pressure head distribution 
below the weir. The effect of joint increases with 
decreasing depth to it from the base of the weir. This 
is due to the fact tnat at shallow depth, the joint 
meets higher equipo tentials and these are shifted 
towards the downstream side. However, it should be 
once again noted that, as it is not possible to 
control the thickness of the stiver paint, simulating 
the secondary permeability absolutely, the results should 
be viewed qualitative than quantitative. 

Figs. 4.4 through 4*7 shov7 the distribution 
of uplift pressure along the base of the weir for 
A = 0.2, 0.4, and 0.6 respectively for various 
values of and , It should be noted that in 1he 
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above figures the values corresponding to j"'! =0.5 
are not plotted as these were lying closely between 
the values presented for = 0 and 1.0. Also in 
Pig. 4.4 the values ofO^ are restricted upto only 5 
as the joint is located at <^= 0.2 and as such the 
minimum value of o' for which results could be obtained 
would be 5.0 so that the sheet pile does not cut the 
joint. Por o( = OQ , the case corresponds to that of 
no sheet pile for which the analysis and results have 
been presented eaT’lier by Madhav and lahshmidhar (1968) 
It is seen from Pigs. 4.4 through 4*7 that witii increa 
sing length of the sheet pile there is considerable 
reduction in the uplift pressure for the pile located 
at the upstream end. However, the trend gets reversed 
when the pile is located at downstream end. 

These results are in confirmity with the already avail- 
able solutions based on the assumption of homogeneous 
porous media. 

Pigs. 4-7 through 4-9 show the effect of 
location of sheet pile on pressure head at D (for 
definition of point D, please refer inset of Pig. 4-7) 
for A. = 0.2, 0.4, and 0.6 and for different values 



42 


of rf . Similar information for point 'C which 
corresponds to the tip of the sheet pile is presented 
in Figs. 4o 10 throu^ 4.11. From symmetry, it can 
be considered that, pressure head at the tip of the 
sheet pile located at the centre of the weir will be 
50 io of head difference irrespective of the length 
of the pile in case of homogeneous porous media. 

Such a trend is seen in Figs. 4.11 and 4.12 where in 
the location of joint is relatively away from the base 
of the weir. However, the pressure head at the tip 
of the sheet pile falls below 50 i° because of the 
presence of a joint in the near vicinity of the weir 
as shown in Fig. 4.10. 

Figs. 4.15 and 4.14 show the effect of position 
of joint on pressure head at h for '\'\= 0 and = 0,5 
respectively for various values of o', . From Fig, 4. 13 
it can be concluded that with increasing distance of 
location of horizontal joint from the base of the weir 
the pressure head at D increases and would eventually 
tend to value corresponding to that of no joint 
case. Also it is seen that for a sheet pile located 
at upstream end the pressure head at D is considerably 
reduced for smaller values of meaning thereby 
a joint at smaller depth, This is an expected 
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behaviour as the pressure would decrease because of 
the presence of draining joint in the vicinity of the 
weir. Hovrever, v?hen Ihe sheet pile is located at the 
centre, the pressure head at D is not markedly influenced 
as can be seen from Pig. 4. 14. 

Pig. 4.15 shows the variation of uplift pressure 

along the base of the weir for the pile ab upstream and 

for various locations of the joint (i.e. for various 

values of A ). This Pig. clearly bring- out the 

necessity of considering the presence of rock jointing 

in the analysis as it i s seen that the pressure 

distribution gets considerably altered. Also as reported 

earlier by Madhav and Lakshmidhar ( 1968) for C< = ch s 

in 

it is seen that there is increase/pressure at downstream 
end because of the presence of a joint in the vicinity 
of the v^eir for all values of -A • 

(ii) v: saTiaAi joint - 

Pigs. 4. 16 and 4.17 show the effect of 
location of joint as influencing the equipo ten tial 
lines for a sheet pile at upstream end of the weir 
for 'V'\ =0.2 and 0.8 respectively. In both these 
Pigs, for the same location of the joint, eq.uipotential 


I 
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lines are shown for various values of CA . As indicated 
in the Pigs » , equipo ten tial lines are considerably 
affected by the presence of joint. 1'lh.en joint is 
located towards the downstream end of the weir the 
effect of joint on equipo ten tial line decreases in 
relation to the joint nearer upstream end. Ey compar- 
ing the equipo ten tial lines for same value of {y , 
it is observed that the location of the joint i.e. 
its relative distance from the upstream end affects the 
pattern of pressure distribution beneath the base of the 
weir.. 

Pigs.. 4. 18 throu^ Pig. 4.20 show the 
distribution of uplift pressure at the base of the 
w.eir for 1'] = 0.2, 0,4 and 0,8 respectively for various 
values of V' ^nd cy. . In the above figures the values 
corresponding to iM = 0.5 are not plotted as these 
'were lying closely between the values presented lor 
i'-i - 0 and 1.0. It is noticed from Pig. 4*17 through 
4.20 that the trend of the influence of cA is same 
as in the case of horizontal joint. 

Pig. 4-21 through 4.23 indicate the effect 
of location of sheet pile on pressure head at- D 
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(for definition of point D, please refer in set of 
Pig. 4*7) for = 0.2, 0^4 and G.8 and for different 
values of :/ (5-0, 2.5 and 1.7). Similar inforraation 
for point ’O' corresponding to the tip of sheet pile 
is presented in Pigs, 4-24 through 4.26. It is seen 
from these fi.gures that the pressure at the tip of the 
sheet pile is not 50 for pile located at the center 
as is true from symmetry for a sheet pile in homogeneous 
porous media. This deviation is becsiuse of the presence 
of joints in the vicinity of the sheet piles. 

Pigs, 4,27 and 4.28 show the effect of location 
of joint on pressure head at D for H = 0 and H^.5 
respectively for various values of yf . Por the pile 
at upstream end it is seen that pressure head at P 
would reach an optimum value for a joint located around 
'Vi 0.5 for sheet piles having C/ = 5 and 2.5, Plowever? 
for the pile located at the center of the weir, the 
minimum pressure head is attained for a joint located 
at 'V'i = 0,5. This is obvious because in this case 
because of pronimity of vertical pervious joint there 
would be considerable reduction in pressure head at D, 

Pig. 4.29 shows the variation of uplift pressure 
along the base of the weir for = 0 and 1.0 for 
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~ 0,2j 0,4 and 0,B, Ihis Fig, clearly brings out 
the necessity of considering hie presence of rock 
jointing in the analysis as it is noticed that the 
pressure distribution gets considerably influenced. 
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TABIiB 4.1(a) HOEIZOMI-ili JOINT 


^ = 0 . 0 , 0 . 5 , 1.0 


A= d/b 

cl- s/b 


.2 

10, 

5 

.4 

10, 

5, 2.5 

. 6 

10, 

5, 2.5, 1.7 


TABLE 4.1(b) 7EETICAL 

JOINT 


U = 0.0, 0.5, 

1.0 



4]- d/b 

a 

-= s/b 


.2 

'5, 

2,5, 1 

.7 

•4 

5, 

2.5, 1 

.7 

,8 

5, 

2.5, 1 

.7 
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CHAPTER 5 


IPPECT OP LOYiPiR PijOV/ BOUHC-YRY AMD. AIJISOTROPY OP TEiE POROUS' 
MEDIA OE THE liPLIPT PRiiiSuLUE POE A lEIR ITITH A SIEGIE 1/ERT- 
lOAE SHEET PIES; 

5.1 GBMBRAE: 

The confined steady flow for a weir with a vertical 

sheet pile resting on homogeneous mass of infinite extent 

has been investigated independently by Khosla et al 

( 1954 ) and Pavlovsky (1933) (ride Harr (1962)). Por the 

same problem treating the porous media to be of finite 
depth and underlain by a horizontal impervious layer, 

solutions have been presented, by Pavlovsky and Muskat 

(vide Ha.rr, 1962). Instead of assuming the underlying 

stratum to be impervious, Grinsky (vide Harr, 1962), 

presented analysis taking the stratum to be a pervious or 

a draining one. In all these cases, the flow domain 

was assumed to be isotropic and homogeneous. However, in 

many instances pervious stratum through which flow is 

taking place is anisotropic in nature. Experimental 

investigations of Graton and Eraser ( 1935), Arnovici 

( 1947 ), Reeve and Kirkhan (1951) Mansur and Dietrick 
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( 1965 ), and Boulton ( 1970) clearly show the maisked 
evidence of anisotropio nature of permeability. Even- 
though, this fact is long- since recognized not much 
work has been done in this area only recently Reddy and 
Misra ( 1972) analysed many problems of confined flow tak- 
ing into account anisotropic nature of soil. However, 
there are many cases for which no solutions exist and as 
such in the present investigation results are presented 
of an electrical analog study dealing with weirs resting 
on anisotropic soil and for various flow boundary 
conditions . 

5.2 STATBlffilT OE DIE PROBESM: 

Problems concerned with anisotropic porous media 
are invariably solved by transforming the actual flow 
domain into a fictitious isotropic region by a suitable 
coordinate transformation, The required scale 
transformation c ai be derived from the equation of 
continuity as follows (Harr, 1962). 


The equation of c© ntinuity for a two dimensional 


steady flow is; 
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ou 

O X 


y 
d y 


0 


(5.1) 


where 

u, T = discharge veloci-fcy in x and y directions, 
respectively. 


Prom the generalised Darcy's law the relationship he tween 
discharge velocity and hydraulic gradient can he written 
as : 


u = 

- (K^) 

- X ^ 

(5.2) 

and V = 

- (%%) 

( ) 
y 

(5.3) 

in Vii'hich: 





K:y = principle coefficient of permeability 
in X and y directions, respectively. 


h = 



p = pressure, 

= unit weight of water and 

x,y = 


co-ordinate axis. 
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Substitution of valuos of u and v in eq.uation 5.1 
leads to; 


K .X 


12 


'h. 


+ K 


-y 


O’ X 


o-y2 


= 0 


(5.4) 


Putting X 


X (^ty/^x) 


eqn. 5.4 reduces 


a \ 


•4“ 



= 0 


(5.5) 


Alternatively, substituting Y = y (^ x/^.y ) "'■^^ 
transforms eq.n. 5-4 into; 


'h 


■■■■■' X 


+ 


d h. 

1 


= 0 


(5.6) 


Thus, by choosing either of the above transformations 
a homogeneous .an isotropic region can be converted into 
a fictitious isotronic region for which the laplace's 
equation is valid. 


Based on above general transformations the follow- 
ing example vide Harr (1962) shows how a weir with a 
vertical sheet pile gets transformed into an inclined 
sheet pile in fictitious isotropic domain. 

1. 1. T. KANPUR ji 

CENim I 

Mi- 2.3 9. 9 6 I 
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Fig. 5.1 represents a cross-section tlirough. a 

row of impervious sheet piles, founded in an anisotropic 

base. The permeability characteristics of the base soil is 

given by an ellipse of permeability inclined at 45° 

to the horizontal with = 4k2. To get the section 

tronsforraed, a coordinate system, parallel to the main 

axis of elipse is established through some point in the 

flow domain (point 2 Fig. 5.1). The boundaries of the 

fictitious flow region, wherein Laplce's equation is valid 

are then obtained by multiplying the perpendicular distance 

from the boundaries of the natural flow domain to the 

reference axes by the approximate scaling factor. In 

Fig, 5.1(a) all the distances parallel to 1-1 ax:is 

were reduced by (k 2 /k-| ) The primed numbers locate 

the vertices of the transformed section. Once the 

solution has been obtained for the transformed section 

(the transformed flow not is given in Fig. 5.1(b), the 

solution for natural boundaries ( Fig .o 1 c ) can be 

obtained by applying the inverse of scaling factor 

(kVk 2 )'^'^^ =2). It is clear from Fig. 5.1(b) that in 

case 

general, for general anisotropic/the vertical sheet pile 
becomes an inclined one in the transformed equivalent 
isotropic region. 
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In addition to the localised non— homogeneity and 
anisotropy of porous media another parameter which influences 
considerably the confined flow beha.viour is the nature of 
flow bound ciry» Becauso of its importance many flow 
boundaries have been considered by various investigators 
in their analysis and Big. 5* 2. summarises the various 
cases for which "theoretical solutions are available. It 
is seen from "this figure that most of the earlier solutions 
were based on the assumption of finite porous media under- 
lain either by impervious or draining stratum. only recently 
Reddy, Mishra and S£efc*¥‘amiah (1972) presented the 
theoretical solutions wherein the impervious flow boundary 
was considered as inclined to the horizontal surface. 

Electrical enalogy me'thod has been also used to 
study .the effect of naturo of lower impervious boundary 
as affecting the uplift pressure. Big. 5*3 summarises "the 
cases studied by this technique by various investigators. 

However, in many field situations the lower flow 

boundary is neither horizontal nor inclined to the 

horizontal stratum at one inclination. One typical field 

cross-section shown in Fig. 5.3j which is quite of 
commo'n occurrence, indicates that it is necessary to 

inves'tigate the case of a weir resting on pervious 
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stratum uiiderla,in by an impervious stratum with 
sloping impervious strata downstream (Pig. 5.4). 

Also this case is most g^^neral in tho sense that if 
the distance m (Pig. 5-4) is very large as compared 
to other dimensions, the case corresponds to that of a 
finite porous media underlain by impervious strata and 

when d is very large as compared to other dimensions, 
the case corresponds to that of a inclined impervious 
boundary. Simulating such a flow boundary results are 
obtained for a weir with one sheet pile resting on 
anisotropic soil. As has been explained .earlier the 
vertical sheet pile becomes an inclined one in the 
transformed section and as such in electrical analogy 
study the sheet pile is taken to be an inclined one to 
account for anisotropy of the soil. It should be noted 
that the results thus obtained for this case and 
subsequent cases are for the transformed isotropic 
zone. 


The second flow boundary which has been consider- 
ed is shown in Pig. 5.5> where it is assumed that 
a sloping draining stratum exists at downstream of the 
weir. Once again tho medium through which flow is 
taking place is assumed to be anisotropic. It should 
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iDG noted here that the .present case falls at the 
other extreme for which Misra et al presented the 
solution. It is obvious that the solution of Misra 

j 4-1 -I j. ^ _^i 

and the results of the present study s.affle as two 
bounds to a two layered soil system seperated by the 
inclined flow boundary. 

The last flow boundary condition considered in 
the analysis is schematically shown in Pig. 5.6. 

Because the underlying and the sloping stratum are 
assumed to be free draining the analysis of this case 
in conjunction with results of 1st case would serve as two 
bounds for a layered soil system separated by the flow 
boundary. 

The different parameter that are varied in 
all the above cases are defined in Pigs. 5*4, 5*5 and 
5.6. Por presenting the results in non-dimensional 
form following dimensionless parameters are defined for 
various cases studied. 

t'- = m/b 

1'* = b/s 
( 

A = d’A 

= b/b 
C/A b/s' 
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where m = distance of the inclined pervious or 

impervious boundary from toe of the weir. 


b = base width of weir, 

s' = embedded length of inclined sheet pile? 

s = embedded length of vertical sheet pile , 

d' = distance of pervious or impervious stratum 
from base of the weir. 


= angle in degree made by embedded length 
of sheet pile with horizontal boundary as 
shown in Pigs. 5.4 and 5.5j 

b^ = distance of sheet pile from heal of the weir. 


The value of .'f will obviously depend upon 
the inclination of tiie axis of major or minor principal 
permeability with horizontal axis. Tables 5.1, 5*2, and 
5.5 summerise. tae values of various v.ariables as 
defined above for which results have been obtained in 
the present study. 
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5*3 ESSUL'iS MD DIBOuSSiul-l uP HiiSULTSs 

Case (i) weir resting on a finite anisotropic pervious 
soil underlain by an impervious stratum with inclined 
sloping impervious boundary? downstreamo 

To study the effect of anisotropy of the soil, 
the angle ) was varied from 30° to 120° with 30 ° 
interval. For studying the -effect of lower flow boundary 
and y were taken equal to 1.0. These values 
have been considered in the analysis in preference to 
many other combinations as the primary aim was to study the 
presence of such a flow boundary in the vicinity of 
location of weir. However the variation in geometry of 
the lower flow boundary is achieved by taking © = 60° 
an d 30°. 

Pigs. 5.7 through 5*14 show^ the variation of 
uplift pressure along the base of the weir for various 
inclinations of sheet pile and downstream 
sloping impervious boundary. Por studying the effect 
of position of sheet pile as affecting the uplift 

distribution two extreme cases namely (i) pile 


pre ssure 
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at upstream end and (ii) pile at downstream end 
have been studied. Fig. 5.7 shows a typical variat- 
ion of uplift pressure at base for 9 = 60°, 'Y = 30° 

A = 1.0 and = 1.0 its can be expected iti is seen 
that with increase in length of the sheet pile there 
is considerable decrease in uplift pressure for a 
sheet pile located at upstream and whereas the trend 
is exactly reverse in case of pile located at the 
downstream end. Similar trends are observed in all the 
Figs. 5.8 through 5.14. The effect of anisotropy of 
soil as affecting the uplift pressure at the base is 
shown in Fig. 5.15(a). It is seen from this figure 
that the uplift pressure distribution is considerably 
influenced by the inclination of major axis of the ellipse 
of anisotropy with horizontal direction. The effect of X ? 
on uplift pressure is shown in Figs. 5.15(b) and 5.15(c) 
for Oi. = 'A- , B = 1.0 and 0 = 60° and 30° respectively. 

For 0 = 60° the results exhibit similar trend as in the 
absence of sloping impervious stratum (Harr 1962). However, 
when 9 = 30 °, implying thereby inclined impervious boundary 
more closer to weir, the trend is changed in relation to 
that for 9 = 60°. 
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Qase Weir resting on infinite anisotropic porous 

mediuni with sloping pervious boundary. 

To study the effect of inclination and the 
location of pervious sloping boundary results have been 
obtained for 9 = 15°, 30° and 60° and = 1.63, 1.0 and 

0.5 for the case of a weir with no sheet pile ( i. e • O'',, = oO) . 
These results are presented in Figs. 5.16 through 5.18. 

Fig. 5.16 shows the variation of the uplift pressure 
distribution along the base of the weir for |h = 1.65 for 
various values of 9. As can be seen from the graph, Ihe 
pressures decrease considerably as 9 decreases. This is 
obvious, beca,use with nearer pervious boundary, the higher 
eq.upo tential lines turn more towards the upstream side of 
the weir thereby 'decreasing the pressures along the 

Cl 

base. Fig. 5-19 shows the effect of 0 and on the 
pressure head at the center of the weir. It is seen from 

i 

this fi.gure that with decreasing k and 9 there is 
considerable reduction in the uplift pressure. 

Fig. 5.20 through 5.27 show the variation of 
uplift pressure distribution along the base of the weir for 
= 30°, 60°, 90°, 120° and 9 = 30° and 60°. These 

results are for =1.0 and for = 5.0, 2.5 and 1.25 

and =0 and 1.0. Fig. 5.20 shows a typical variation 
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of uplift pressure distribution along the base of the 
weir for 9 30°, 30°, B= 1 «0, =0 and 1.0 

and '.A = 1 . 25 } 2 . 5 ? and 5-0. Ihe figure once again 

brigs out the effectiveness of placing the sheet pile at 
the upstream end for reducing the uplift pressures. 

Pig. 5 .28 and iig. 5.29 show the effect of location of 
sheet pile on pressure head at tip of the pile for m/b =1.0 
and 0.5 respectively ior both © = 60 ° and 30° and for'f 
equal to 90°. These results show that with decreasing 
value of m/b there is considerable decrease in pressure 
head at the tip of the pile. Comparison of curves for same 
m/b ratio and tiut with 0 = 30° and 60 ° shoves that 

wm th decreasing value of angle © there is considerable 
decrease in pressure heads. Similar results for point 
D (Refer inset diagram) are presented in Pig.5,50and 5.31 
exhibiting similar trends. 

C ase (iii) weir resting on finire pervious anisotropic soil 
underlain by draining stratum and whth inclined pervious 
boundary downstream; 

In the first instance, to study the 
effect proximity of such a pervious, boundary as affecting 
the uplift pressure, results were obtained lor/\= 1.0 and 
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0.5 for B ^ 1.0, as 0 - 60° ard 30°. Pigs. 5,32. 

and 5033 show ih© corresponding results.Ii ig seen from 
these figures that with decreasing values of X 
(thereby implying nearer draining stratum ) and Q 
the pressures decrease considerably along the base of 
the weir. To study the combined effect of anisotropy 
and lower 'flow boundary results have been obtaindd 
for = 30 °, 60°, 90 ° and 120°, = 1.0 and X= 1.0. 

The variation in lower f low boundary being achieved 
by varying angle 0 as 30° and 60°. These results 
are presented in I'igs. 5 <34 through 5.41. I’ig. 5.34 
shows a typical variation of uplift pressure along the 
base of the weir for © = 60° and {= 30 * for various 
values of (/ and for pile located either at upstream 
or downstream end. By comparing Pig. 5.34 of this 
case to Pig. 5»7 of case (i) it is seen that because 
of the extreme change in the nature of the flow 
boundary from impervious to a pervious one there is 
considerable decrease in uplift pressure distribution. 
>Similar trends can be also observed by comparing 
Pigs. 5.8 through 5.14 with 5.35 through 5*41. 
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T£BIB 5.1 IMPERVIOUS SIPuATM TittlH SIOPIUG- 
IMPERVIOUS STRATUM. 

V' = 0, 1.0 


(=m/b) x\'(d'/b) y a', or x' 





O 

o 

= 5,2.5,1.25 




60° 

= 5,2.5,1.25 




90° 

= 5,2.5, 1.7 



60° 

120 ° 

= 5,2.5, 1.25 

1.0 

1.0 







O 

O 

= 5,2.5, 1.25 




60° 

= 5,2.5,1.25- 



50° 

O 

o 

= 5,2.5,1.70 




120 ° 

= 5,2.5,1.25 


TABLE 5.2 PERVIOUS SLOPING BOUNDARY 
5.2(a) Nc Sheet Pile 


0 

^ ( -rm/b ) 


.5 

60° 

1.0 


1.65 


.5 

50° 

1.0 


1.65 


.5 

15° 

1.0 


1.6 
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TA13EE 5.2 (b) TEHTICAL SH3:3T PIIS 


f-i ( - b /b ) 

= 0, 0.5, 

1.0 

e 

t5 

On 

o 

O 

1.0 

5,2.5, 1.7 


0. 5 

5,2.5, 1.7 


1.0 

5,2.5, 1.7 

O 

O 

0.5 

5,2.5, 1.7 

LABIS 5.2(c) 

iiciriiroD 

SHBET Pir® 


= 0 , 1.0 B = 1.0 


G 

■ 7 : 



30 ° 

5 , 2 . 5 , 1.25 

o 

o 

0 

0 

VO 

5 , 2 . 5 , 1.25 


120^ 

5 , 2 . 5 , 1.25 

°o 

30 '^ 

5,2.5, 1 .25 

0 

0 

5 , 2 . 5 . 1.25 


0 

0 

C\i 

5 , 2 . 5 , 1.25 
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TASIE 5.. 3 PERVIOUS ST'RAIUM WITE BTOLIESD PERVIOUS 

BOUITOARY 

\-\ = 0 , 1,0 

. = 1.0 

X i o 


, 'y O' 5 ' 


30 ° 

5 , 2 , 5 , 1.25 

o 

o 

60 ° 

5 , 2 . 5 , 1.25 


90 ° 

10 , 5 , 2.5 


120 ° 

5 , 2 . 5 , 1,25 


30 ° 

5 , 2 . 5 , 1.25 

O 

O 

60 ° 

90 °. 

5 , 2 . 5 , 1.25 


10 , 5 , 2.5 


120 ° 

5 , 2 . 5 , 1.25 
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G-^NSPAI CaiGIUSIONS 


Based on the studies reported in Chapters 4 and 5 , 
the following broad based conclusions are arrived at? 

(i) It is seen that by adopting few modifications to 
the standard electrical analogy methods the problem of 
conf filed flow through rock mass can easily be solved. The 
secondary permeability due to flow of water being easily 
simulated by dravdng a line of silver paint on teledeltos 
paper. 

(ii) The results indicate that the uplift pressure is 
influenced considerably by the presence of either vertical 
or horizontal joints. The degree of influence of joint 
being dependent on the disposition of joint in relation 

to the base of the weir and location of sheet pile. 

(iii) The pressure head immediately next to sheet pile 
located at upstream end is considerably reduced when the 
horizontal joint is located near the base of the weir. 
However, there is not such a marked influence of joint on 
pressure distrihution immediately next to sneet pile wnen 
the pile is located at the center of the weir (Pig. 4.14 and 

4.15). 

(iv) The pressure head at the tip of the sheet pile 
is reduced because of the major horizontal joint near 
base of the weir (Pig. 4.10). 
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(t) The pressure head immediately next to sheet 
pile is seen to be considerably dependent upon the 
relative position of sheet pile and the joint along the base 
of the weir (Figs. 4.27 anci 4.28). 

(vi) It is seen that the anisotropy of the soil as well 
as the nature of the lower flow boundary have considerable 
influence on the uplift pressure distribution (results of 
Chapter 6). 

(vii) The results obtained for various lower flow bound arj?' 
conditions can serve as lower and upper bounds for a two~lay^€tA. 
soil system separated by a boundary corresponding to the 
geometry of the lower flow bound .ary. 

Because of the versatid.ity of the electrical 
analogyme thod problems of confined and unconfined flow 
through the jointed rock mass for many practical situations 
can be analysed. Also the present study can be extended 
for the case of confined flow belovj weir with multiple 
sheet pile resting on anisotropic soil. 
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